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The density perturbations generated when the inflaton decay rate is perturbed by a light scalar 
field x are studied. By explicitly solving the perturbation equations for the system of two scalar 
fields and radiation, we show that even in low energy-scale inflation nearly scale-invariant spectra 
of scalar perturbations with an amplitude set by observations are obtained through the conversion 
of x fluctuations into adiabatic density perturbations. We demonstrate that the spectra depend 
on the average decay rate of the inflaton & on the inflaton fluctuations. We then apply this new 
mechanism to string cosmologies & generalized Einstein theories and discuss the conditions under 
which scale-invariant spectra are possible. 
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In inflationary cosmology the seeds of large-scale struc- 
ture are originated from quantum fluctuations of a light 
scalar field responsible for inflation- called the infla- 
ton. The adiabatic perturbations generated from inflaton 
fluctuations exhibit nearly scale-invariant spectra, as re- 
quired by observations [1] . The energy scale of inflation is 
determined by the amplitude of large-scale perturbations 
in the Cosmic Microwave Background (CMB). For exam- 
ple, in the simple model of chaotic inflation with potential 
V = (l/2)m^(j> 2 , the mass of the inflaton <f> is constrained 
to be ~ 10 _6 m p i by the Cosmic Background Explorer 
(COBE) normalization [2,3] (Here m pl ~ 1.22 X 10 19 GeV 
is the Planck mass). 

On the other hand, if inflation occurs at a much lower 
scale, the amplitude of density perturbations is too low to 
explain the temperature anisotropy in CMB. In this case 
we need to consider some other mechanism to generate 
sufficient density perturbations. 

Recently a number of authors [4,5] independently pro- 
posed an alternative method to produce scalar metric 
perturbations using the spatial variation of the inflaton 
decay rate, T. Since the coupling strength is generally a 
function of the vacuum expectation value of scalar fields 
in string theory, it is natural to consider the case where 
the coupling fluctuates due to scalar field perturbations. 
When there exists some light field whose mass is smaller 
than the Hubble rate (such as modulus [5]), its large- 
scale fluctuations are not exponentially suppressed dur- 
ing inflation. This can provide another source of scalar 
metric perturbations in addition to the inflaton. In fact 
the authors in [4] considered the decay of the inflaton 
by treating it as a matter fluid and showed that density 
perturbations can be generated by the perturbed decay 
rate, ST. 

In this work we shall extend the analysis of Ref. [4] to 
the more realistic situation where the inflaton is treated 
as a scalar field, and completely follow the evolution of 
metric perturbations during inflation and reheating. This 
is particularly important when we treat the dynamics of 
reheating consistently, since the fluid approach can lead 



to some loss of information (like parametric resonance) 
by averaging over the scalar field oscillations. We will an- 
alyze the case where perturbations of both the inflaton 
and the decay rate coexist, and estimate the amplitude 
and the spectral index of metric perturbations after in- 
flaton decay. We shall also apply this new mechanism to 
string-inspired cosmologies & generalized Einstein theo- 
ries. 

Let us consider the chaotic inflationary scenario with 
potential V = (1 / 2)m 2 ^(j> 2 . During reheating, the inflaton 
(j> decays to standard light particles through the interac- 
tion, Ao^ccr, whose decay rate, T, is proportional to Aq 
[2]. The inflaton can also decay via non-renormalizable 
interactions with superfields, e.g., <j){q/M)qq (M is some 
mass scale), in which case the vacuum expectation value 
of a scalar component of the q superfields, (x), gives rise 
to an effective coupling, A = (x)/M [4]. Following Ref . [4] 
we shall analyze the case where the effective coupling 
A(x) is dependent on the field x as 

X(x)=Xo(l + x/M + ---) . (1) 

When the coupling, (l/2)g 2 (j> 2 x 2 , exists during the re- 
heating stage, this can lead to explosive particle pro- 
duction called preheating [6,7]. In order for preheat- 
ing to occur, we require the large resonance parameter, 
q IC = g 2 <j) 2 /(Am^) ^> 1, at the end of inflation [7], in 
which case large-scale perturbations in x are exponen- 
tially suppressed during inflation due to the heavy effec- 
tive mass, g\4>\, relative to the Hubble rate, H . Then the 
perturbations of the coupling A(x) are vanishing small, 
which do not affect the evolution of metric perturba- 
tions. Therefore in this work we do not take into ac- 
count the effect of preheating coming from the interac- 
tion, (l/2)g 2 (j> 2 x 2 ■ When this interaction is added to the 
Born decay term, Ao<^x 2 j this gives rise to the form of the 
effective coupling (1) perturbed by the inflaton <f> instead 
of x- However we are interested in the case where the 
coupling is perturbed by a light field x other than the 
inflaton. 

Hereafter we will analyze the two-field system of 4> and 
X with radiation. We assume that the field x nas a small 
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mass, m x , relative to the Hubble rate, H = a/ a, during 
inflation (here a is a scale factor). Radiation is generated 
through the decay of the inflaton with the x-dependent 
coupling (1). Consider the following perturbed metric for 
scalar perturbations in the longitudinal gauge [9]: 



ds 2 = -(1 + 2$)dt 2 + a 2 (t)(l - 2^)5 ij dx i dx j 



(2) 



At linear order one has $ = due to the vanishing 
anisotropic stress. Then the Fourier transformed, lin- 
earized Einstein equations are [10]* 

(fc 2 /a 2 )$ + 3H(<f> + H<f>) - (47t/to 2 1 )(0 2 + x 2 )$ 



= - (4tt/ m 2 !) <£5<fi + m 2 ^(pS(f) + xSx + m x xSx + S/-J. 
54> + {3H + T)5<p + (k 2 /a 2 + m£) 5<p 
= - 2m!^$ - T0(2$ + 5T/T), 



5 X + 3HSx+(k 2 /a 2 



Sx = 4x < l' - 2m 2 x4> 



5(i + AHSfj, = 4/x $-(fc/3a)u 



250 + (<5r/r)0 



ii + (4i? + (i/n)v = (k/a) <£ + (5n + 3r^)/4/z 



(3) 

(4) 
(5) 
(6) 

(7) 



where k is a comoving wave number (the k subscript is 
implicit) and T is the inflaton decay rate that is smaller 
than the Hubble rate during inflation. 5/j, is the pertur- 
bation of radiation whose energy density is given by fi, 
and v is the velocity of the radiation fluid. 

The decay rate, T, is proportional to A 2 , and A is de- 
pendent on the field x through the relation (1). Therefore 
the fluctuation of the decay rate can be written as 



<fi can be the same order during inflation [15]. This cor- 
responds to double inflationary scenarios driven by two 
scalar fields, in which case the strength of the correlation 
is numerically investigated in Refs. [15,16] for 5T = 0. 

In the presence of the 5T fluctuation, the situation is 
different from the one discussed in Ref. [14]. First of all, 
even if \ is much smaller than 4>, the conversion of the 
ST fluctuation into the gravitational potential occurs as 
seen from Eqs. (3), (4) and (6). As long as the mass of x 
is light and large-scale x fluctuations are not suppressed 
during inflation, this naturally leads to the conversion to 
the gravitational potential through the perturbed decay 
rate, ST. The advantage of this scenario is that we do not 
need to rely on two x-dependent source terms in Eq. (3) 
in order to obtain sufficient isocurvature perturbations. 
In this sense it is also different from metric preheating 
[17] where the growth of the x-dependent terms sources 
the gravitational potential. 

In the absence of the <5r term, the amplitude of 
scalar perturbations gener ated from inflaton fluctuations 
are estimated as Sh = \/5127r/75 V 3/2 / (ra^l \ ) \ k= a H 
[3]. Substituting the value <p — 3to p i corresponding 
to the 60 e-foldings before the end of inflation, we get 
Sr — I2m^/m v \. Making use of the COBE result 
Sr — 2 x 10~ 5 , the inflaton mass is constrained to be 
77i0 ~ 1.7 x 10~ 6 m p i. The spectral index n for the power 

spectrum, V<s> = ^-(\^k\ 2 ), is estimated as [3]: 



1 - 6e + 2t, ~ 1 
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~ ^1* and 77 = r H^.Y!L _ 

are the slow-roll parameters. This is a slightly red-tilted 
spectrum for -C H. 

The above basic picture of the generation of density 
perturbations can be substantially modified when the ST 
fluctuations are involved. Let us first estimate the spec- 
tral index of the S\ fluctuation in Eq. (5) with mass satis- 
fying 77i x < H. When the field \ is subdominant relative 
to <f> (x <§C 0), we can neglect the r.h.s. of Eq. (5) coming 
from the gravitational potential. Then the spectral index 
for the super-Hubble 5\ fluctuation is estimated as [18] 



where we dropped the contributions coming from the 
higher-order terms in Eq. (1). When the mass of the 
field x is smaller than iJ, large-scale % fluctuations are 
not suppressed during inflation, thereby providing source 
terms for the gravitational potential $. 

We shall analyze the situation where the background 
field x is not dynamically important relative to the in- 
flaton <j) in order to avoid the second stage of inflation 
to occur. This means that x "C </> in Eq. (3), in which 
case the correlation between adiabatic and isocurvature 
perturbations is small [14] due to the negligible angular 
velocity in field space (tan 5 = x/<j) — 0). In the absence 
of the ST term the conversion from isocurvature to adia- 
batic perturbations occurs significantly only when x and 



3H- 



■(m 2 - to 2 ). 



(10) 



*After this work was completed, the authors in Ref. [11] de- 
veloped a gauge-invariant formalism for the large-scale curva- 
ture perturbations with spatial and time variations of the in- 
flaton decay rate using the formalism in Ref. [12]. See Ref. [13] 
for a related work. 



In the presence of the ST fluctuation, the power spec- 
tra V<s> exhibit a spectral tilt which is a mixture of 
Eqs. (9) and (10). When m x is smaller than to^, one has 
slightly red-tilted spectra of $ which are close to scale- 
invariant (see Fig. 1). The spectra of Sx are blue-tilted 
for m x > rritp, but this does not necessarily mean that the 
gravitational potential is also blue-tilted. In fact, even 
when to x > 3rn0, super-Hubble Sx fluctuations begin to 
be exponentially suppressed during inflation, in which 
case the gravitational potential is dominated by the in- 
flaton fluctuation whose spectrum is given by Eq. (9). 
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Typically we have red-tilted, nearly scale-invariant spec- 
tra of $ for wide range of the parameter space. 

If the inflaton mass is less than ~ 10 _6 m p i, the am- 
plitude of the power spectra V<s> is suppressed compared 
to the observed value, ~ 10~ 9 (see Fig. 1). Taking 
into account the ST fluctuation, the amplitude of density 
perturbations can be much higher. In fact, as shown in 
Fig. 1, it is possible to get nearly scale- invariant power 
spectra whose amplitude are of order 7^$ ~ 10 -9 . 
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FIG. 1. The spectra of the gravitational potential $ af- 



ter inflaton decay for m</> = 10 m p i, 



= 0.87710, 



O.lm^ with initial conditions 

6„ 



M = 1.2 x 10"°ro P i, r 

4>i — 3.2m p i and x> = 10 _t, m p i. Nk is the number of 
e-foldings before the end of inflation. The solid and dotted 
curves correspond to the case with ST and 8T — 0, respec- 
tively. We find that the required power spectra, V<s> ~ 10 -9 , 
are obtained by taking into account the ST fluctuation. In- 
set: The evolution of Pj> for <5r ^ for the mode that crossed 
the Hubble radius around Nk = 60. 
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FIG. 2. The parameter space where the power spectra, 
Vg> — 10 -9 , are obtained for the mode corresponding to 
Nk = 60 in the case of m x = O.lm^. The condition, \ ^ M, 
is satisfied in all cases of this plot. Each curve corresponds 
to the case of the average decay rate, f = T/m^ = 0.1, 0.01, 
and 0.001, respectively. 

We also plot in Fig. 2 the parameter range with re- 
spect to ro^/mpi and ST/T in which the power-spectra, 



V$> ~ 10~ 9 , are obtained after inflaton decay. When 
the x mass is light relative to the Hubble rate, one gets 
the perturbations of order 5x ~ H after Hubble radius 
crossing. Therefore the fluctuation of the decay rate is 
approximately estimated as ST/T ~ H/M for \ ^ M 
from Eq. (8). We can have a considerable amount of fluc- 
tuations in ST/T by varying the mass scale, M. In Fig. 2 
we find that inclusion of the perturbed decay rate allows 
for the wide possibility to generate the perturbations of 
order V<$> ~ 10 -9 even for <C 10~ 6 m p i. This figure 
corresponds to the one where x is much smaller than <f> 
(and M). We wish to stress that 8x fluctuations can be 
sufficiently transferred to adiabatic density perturbations 
even for tan 9 ~ 0. When the average decay rate T gets 
smaller, we require larger values of ST/T (corresponding 
to smaller values of M) to lead to sufficient conversion as 
confirmed by Eqs. (4) and (6). This behavior is actually 
seen in our numerical simulations of Fig. 2. In the limit 
of T — > (Ao — > 0), the above conversion does not occur 
apart from the standard one discussed in Ref. [14]. 

We can apply the above new mechanism for the more 
general action 



S = 



^i?-i(V^) 2 -FM(V X ) 2 -C/ 



(11) 



where F(ip) is a function of ip and U is a potential of 
scalar fields including their interactions. The action (11) 
contains not only various generalized Einstein theories 
[19,20] but also the tree-level Pre-Big-Bang (PBB) cos- 
mology [21,22] from the dimensional reduction and 4D 
effective description of ekpyrotic cosmology [23,24]. 

The effective potential U typically includes an expo- 
nential term in generalized Einstein theories, which can 
lead to inflationary solutions [19,20]. Let us consider 
the soft inflation model [19] with U = Uq exp(— ^Ipptp) 
and F = (l/2)e° v with j3 being a constant (hereafter 
we set the units such that Sn/rn^ = 1). We shall as- 
sume that the decay rate of the inflaton is perturbed 
by the field x with a negligible mass. The evolution of 
the background is given as a a t p oc (-rf^-P) and 
e 1 ? cx {-(p- i)t}v^p/(i-p) w i t h T = J a - x dt{< 0) be- 
ing a conformal time [25] . Assuming that the background 
field x is not dynamically important, the perturbed equa- 
tion for §x is approximately written as 



Sx + (3H + ptp)8x + (k 2 /a 2 )S X - . 



(12) 



Introducing new variables, b = e^ v ' 2 a and Sx = bSx, we 
get 



Sx" + {k 2 - b"/b) S X = 0, 



(13) 



where a prime denotes the derivative with respect to r. 
When the evolution of b is characterized by 6 cx (— r) 7 , 
the spectral index of the Sx fluctuation is given as ng x — 
4 — |1 — 27 1 [26]. In the present case 7 is given by 
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•"—fM'WD- (14 ' 

In the limit of p — > oo , we have 7 — > — 1 , thereby yielding 
the scale- invariant spectra, n$ x = 1. In the presence of 
the noncanonical kinetic term (0^0), we have n,5 x = 1 
for /3 = — \/2/p or (3 = — (3p — 2)y / 2/p (corresponding 
to 7 = —1 or 7 = 2). If the decay rate of the inflaton is 
perturbed by the field x> it is possible to generate scale- 
invariant spectra of density perturbations in a similar 
way as discussed previously. 

In PBB and ekpyrotic cosmologies, one has the con- 
tracting phase of the universe in the Einstein frame before 
the bounce [27]. Let us consider the ekpyrotic scenario 
whose potential is given by U = —Uq exp(—w2/p~(p) 
with Uo < and < p < 1. In this case the 
contracting phase is characterized by the scale factor, 
a oc (-t) p oc (-t)p/< 1 -p), with t < 0. Since the brane 
modulus, tp, evolves as e v oc {(1 — p)(— r)} v ^/( 1_p ), we 
get the same 7 given by Eq. (14) for the fluctuation 
6\- In the case of minimal coupling (0 = 0), one has 
ng x = (3 — p)/(l — p) — 3 for p ~ 0. This is the same 
blue spectral index as the curvature perturbation, 1Z, ob- 
tained in the system of the brane- modulus only [28] . We 
can have ns x = 1 for p = 2/3 [29], but the system is 
found to be unstable in this case [27] . When = — \j2jp 
or — (2 — 3p) y/2/p, metric perturbations can be scale- 
invariant provided that 5\ fluctuations are sufficiently 
converted to <& through the perturbed decay rate. 

The PBB scenario corresponds to the vanishing po- 
tential with p = 1/3. It is convenient to use the re- 
lation <p = ~V2<P between the dilaton <fi in the PBB 
case and the field tp in the ekpyrotic case [27]. This is 
equivalent to replacing (3 in Eq. (11) for new beta 0, as 
(3 = —y/2/3. Then the spectral index for 5\ is given by 
ns x = 4 — \/3|/3|, which allows for scale-invariant spec- 
tra for \0\ = x/3. Since 1/31 = 2 for the standard axionic 
coupling, we have n$ x ~ 0.54 in this case [22]. If we 
take into account the contribution of the modulus ki- 
netic term, -n(Va) 2 , in the action (11) (here n is the 
number of compactificd dimensions), it is possible to ob- 
tain scale-invariant spectra even for = 2 [22]. When 
the modulus a evolves as e a oc (— t) s , we get the spec- 
tral index ng x — 4 — \0\V3 — 2ns 2 . The scale-invariant 
spectra follow for \0\y/3 — 2ns 2 = 3. 

Of course one must make a numerical analysis in order 
to check the efficiency of the conversion to adiabatic den- 
sity perturbations. In string-inspired cosmologies with a 
contracting phase, noncanonical kinetic terms (0 7^ 0) 
are essentially important to yield scale-invariant spec- 
tra in <f> through the perturbed decay rate coming from 
S\- The correlated adiabatic and isocurvature pertur- 
bations were discussed in Ref. [25] for the same action 
(11) when ST = 0, implying the possibility to generate 
scale-invariant spectra in entropy field perturbations for 
< p < 1. In addition to this, the perturbed decay 



coupling provides us a new source for isocurvature per- 
turbations, whose conversion to adiabatic perturbations 
is efficient even when the background field \ 1S dynami- 
cally negligible. 

In PBB and ekpyrotic cosmologies, we have radiation 
just after the bounce through the decay of scalar fields. 
The spectra of perturbations can be affected by the de- 
tails of the background evolution around the bounce as 
was recently analyzed in Ref. [30]. It is certainly of in- 
terest to investigate the final spectra of $ for <5r 7^ 
by passing through the bounce numerically. Not only in 
the case of standard inflation, this can also provide us 
an exciting possibility to generate nearly scale-invariant 
density perturbations with the amplitude required by ob- 
servations in the context of string cosmologies. 
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